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Abstract 

We  have  studied  thermal  expansion  of  the  surface  layers  of  the  hexagonally  reconstructed  Au  (001) 
surface  using  a  classical  Molecular  Dynamics  (MD)  simulation  technique  with  an  Embedded  Atomic 
Method  (EAM)  type  many-body  potential.  We  find  that  the  top-most  hexagonal  layer  contracts  as 
temperature  increases,  whereas  the  second  layer  expands  or  contracts  depending  on  the  system  size.  The 
magnitude  of  expansion  coefficient  of  the  top  layer  is  much  larger  than  that  of  the  other  layers.  The 
calculated  thermal  expansion  coefficients  of  the  top-most  layer  are  about  -4.  93  '  10-5  A/K  for  the  (  262 
'  227  )  A  cluster  and  -3.05  1 0  -5  A  /  K  for  (101  '  87)  A  cluster.  The  Fast  Fourier  Transform  (FFT)  image 
of  the  atomic  density  shows  that  there  exists  a  rotated  domain  of  the  top-most  hexagonal  cluster  with 
rotation  angle  close  to  1°  at  temperature  T  <  1000K.  As  the  temperature  increases  this  domain  undergoes  a 
surface  orientational  phase  transition.  These  predictions  are  in  good  agreement  with  previous 
phenomenological  theories  and  experimental  studies. 
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THERMODYNAMIC  BEHAVIOR  OF  NANO-SIZED  GOLD  CLUSTERS  ON  THE  (001)  SURFACE 

Sun  M.  Paik1,  Sung  M.  Yoo2,  Min  Namkung3,  and  Russell  A.  Wincheski3 

Abstract.  We  have  studied  thermal  expansion  of  the  surface  layers  of  the  hexagonally  reconstructed  Au 
(001)  surface  using  a  classical  Molecular  Dynamics  (MD)  simulation  technique  with  an  Embedded  Atomic  Method 
(EAM)  type  many-body  potential.  We  find  that  the  top-most  hexagonal  layer  conUacts  as  temperature  increases, 
whereas  the  second  layer  expands  or  contracts  depending  on  the  system  size.  The  magnitude  of  expansion 
coefficient  of  the  top  layer  is  much  larger  than  that  of  the  other  layers.  The  calculated  thermal  expansion  coefficients 
of  the  top-most  layer  are  about  -4.93  x  10  5  A  /  K  for  the  ( 262  x 227 )  A  cluster  and  -3.05  x  10  3  A  /  K  for 
(101  X87)  A  cluster.  The  Fast  Fourier  Transform  (FFT)  image  of  the  atomic  density  shows  that  there  exists  a  rotated 
domain  of  the  top-most  hexagonal  cluster  with  rotation  angle  close  to  1°  at  temperature  T  <  1000K.  As  the 
temperature  increases  this  domain  undergoes  a  surface  orientational  phase  transition.  These  predictions  are  in  good 
agreement  with  previous  phenomenological  theories  and  experimental  studies. 

Key  words.  Molecular  Dynamics  (MD)  simulation,  surface  reconsUuction,  surface  orientational  phase 
transition,  surface  thermal  expansion 

Subject  classification.  Computational  Nano-Materials 

1.  Introduction.  In  last  few  decades  there  have  been  extensive  studies  of  surface  reconstruction  and  phase 
transitions  on  noble  metals  and  semiconductor  surfaces  both  theoretically  [1-3]  and  experimentally  [4-6].  It  is  now 
well  known  that  the  top-most  layer  of  the  Au  (001)  surface  is  reconsUucted  into  a  hexagonal  lattice  with  (5 xm) 
reconsUuction  patterns  where  m  is  an  integer  between  20  and  30.  The  reconstructed  layer  is  slightly  (-3.5%) 
contracted,  distorted  and  buckled.  The  rotated  domains  with  rotation  angles  9  =  1°  are  found  experimentally  [7-9]. 
These  rotated  domains  undergo  orientational  phase  transitions  at  the  temperature  around  T  ~  0.757',,,  where  T,„  is  the 
bulk  melting  temperature. 

In  this  paper,  we  use  the  MD  simulation  technique  to  study  the  surface  properties  and  calculate  the  thermal 
expansion  coefficients  of  the  Au  (001)  surface.  We  find  that  the  interlayer  distance  between  the  top-most  hexagonal 
layer  and  the  second  layer,  dj,  decreases  as  the  temperature  increases,  whereas  between  the  second  layer  and  the 
third  layer,  d2.  decreases  or  increases  depending  on  the  size  of  the  cluster.  The  calculated  thermal  expansion 
coefficients  for  the  interlayer  distance  are  about  -3.05  x  10~5  a  /  K  for  (101  X87)  A  cluster  and  -4.93  x  10~5  AIK 
for  (262  X227)A  cluster.  The  Fast  Fourier  Transform  (FFT)  image  of  the  atomic  density  shows  that  there  exists  a 
rotated  domain  of  the  top-most  layer  with  rotation  angle  close  to  1°  at  the  temperature  T  <  1000K.  For  the 
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temperature  above  1000K  the  unrotated  domain  is  found  to  be  stable.  These  findings  are  in  good  agreement  with 
previous  theoretical  studies  [10-13]  and  experimental  results  [7-9]. 

2.  Method.  We  use  a  classical  Molecular  Dynamics  (MD)  simulation  technique  with  a  modified  Verlet 
algorithm  [14]  to  simulate  surface  structures  and  dynamics  as  a  function  of  temperature.  The  Embedded  Atomic 
Method  (EAM)  type  many-body  potential  [15-17], 

£  =  E^[P,]  +  iE'/(r(j).  <2» 

i  ^  i*j 


is  employed  for  interatomic  interactions  between  Au  atoms  where  the  /*’[p]  is  the  embedding  function,  V( rt])  is  a 


two-body  interaction  potential,  and 


The  charge  density. 


P,-  =  Ep>§).  (2-2) 

m ) 

is  the  sum  of  all  atomic  charges  from  neighboring  atoms  j  at  the  position  of  ith  atom.  The  spline  fitted  functions  and 
parameters  for  Au  atoms  can  be  found  elsewhere  [17]. 

Our  system  consists  of  one  fixed  (001)  layer  and  two  movable  (001)  layers.  On  top  of  these  layers,  one 
hexagonal  movable  layer  with  1225  ~  8281  atoms  is  placed.  We  choose  the  hexagonal  top-most  layer  because  our 
earlier  study  [10-11]  showed  that  the  hexagonal  top-most  layer  is  the  most  stable  one  in  Au  (001)  surface.  The  sizes 
of  the  second  layer  and  below  are  much  larger  than  the  top-most  layer  so  that  the  boundary  of  these  layers  does  not 
affect  the  top  layer  atoms.  The  usual  periodic  boundary  condition  is  not  applied  in  this  simulation  to  avoid  unwanted 
boundary  effects  [10-1 1]. 4 

For  a  given  temperature  the  lattice  constant  of  the  fixed  layer  is  set  to  the  bulk  lattice  constant  at  the 
corresponding  temperature,  and  the  movable  layers  are  equilibrated  to  the  required  temperature  for  about  10,000 
time  steps.  After  the  equilibrium  process  the  thermodynamic  quantities  are  calculated  for  another  10,000  time  steps. 

For  the  sake  of  convenience  we  use  the  MD  units  in  the  simulation  where  e=3.78eV,  c=2.89A,  and 
to=(ma2/£)]n  are  chosen,  respectively,  as  the  units  of  the  energy,  length  and  time.  Here  m  is  the  mass  of  a  gold  atom. 
The  one  MD  time  step  is  given  by  A  t  =0.001  to  ,  where  t0  =2.13xl013sec  for  Au  atoms. 


3.  Results.  In  Fig.  1  we  plot  the  interlayer  distance  between  the  hexagonal  layer  and  second  square  lattice 
layer,  dj,  as  a  function  of  temperature  for  the  cluster  size  of  (262  x 227 )  A  (solid  circles)  and  (101  X87)  A  (open 
circles).  The  layer  distances  for  both  clusters  decrease  as  the  temperature  increases.  The  calculated  thermal 
expansion  coefficients  in  the  direction  normal  to  the  surface  are  about  —3.05  x  1 0  A  /  K  for  the  smaller  cluster  and 
-4.93  x  10  5  AIK  for  the  larger  cluster. 

Fig.  2  shows  the  interlayer  distance  between  the  second  (square  lattice)  layer  and  the  third  layer  as  a 
function  of  temperature  for  the  cluster  size  of  (262  X227 )A  (solid  circles)  and  (101  X87)  A  (open  circles).  The 


4  The  periodic  boundary  may  not  be  applicable  to  this  system  at  all  since  our  system  has  two  different 
symmetries  (hexagonal  symmetry  on  the  top  layer  and  square  symmetry  on  the  second  layer  and  below). 
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layer  distance  d2  for  the  smaller  cluster  (open  circles)  increases  slightly  while  it  decreases  for  the  larger  cluster  (solid 
circles)  as  the  temperature  increases.  However  the  overall  changes  are  much  smaller  than  that  of  the  d /.  Here  we 
note  that,  in  actual  simulation,  the  second  layer  (and  below)  are  much  larger  than  the  top-most  layer  cluster  to  avoid 
edge  effects,  but  the  d2  is  calculated  by  taking  the  average  of  the  z  distance  over  the  atoms  that  are  beneath  the 
hexagonal  layer.  The  rest  of  the  second  layer  (i.e.,  the  part  in  the  second  layer  without  the  hexagonal  layer  on  top)  is 
contracted  more  than  that  for  atoms  with  the  hexagonal  layer  on  the  top. 


Figure  1:  The  interlayer  distance  d2  as  a  function  of 
temperature  for  the  cluster  size  of  ( 262  X  227 )  A  (solid 
circles)  and  (101X87)  A  (open  circles).  The  solid  and 
dashed  lines  are  the  linear  least  square  fits. 


Figure  2:  The  interlayer  distance  d2  as  a  function  of 
temperature  for  the  cluster  size  of  ( 262  x  227 )  A  (solid 
circles)  and  (101  X87)  A  (open  circles). 
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Figure  3:  The  FFT  image  of  atomic  density  for  the 
top-most  layer  cluster  of  size  (101X87  )  A  at 
temperature  T=300K.  The  dotted  lines  are  the  guides 
to  the  eye  for  rotation  angle.  The  peaks  in  the  white 
circles  represent  the  hexagonal  peaks. 


Figure  4:  The  density  plot  of  FFT  image  of  atomic 
density  for  the  cluster  size  of  (101X87)  A  at 
temperature  T=1200K. 


In  Fig.  3  we  show  the  density  plot  of  the  Fast  Fourier  Transform  (FFT)  of  the  atomic  density  for  the  top¬ 
most  layer  cluster  of  (101  X87)  A  at  temperature  T=300K.  The  hexagonal  symmetry  peaks  are  clearly  shown  in  the 
bottom  half  of  the  figure.  Because  of  the  square  lattice  underneath  the  cluster  the  square  symmetry  peaks  are  also 
visible  in  the  figure.  We  notice  that  the  hexagonal  peaks  split  in  two  peaks  (see,  for  example,  inside  of  the  white 
circles  in  the  figure)  indicating  that  the  rotated  domains  coexist  with  unrotated  domains.  The  separation  of  the  two 
peaks  is  very  small,  i.e.,  the  rotation  angle  is  very  small.  The  rotation  angle  is  about  9=1°.  These  rotated  domains 
may  undergo  an  orientational  phase  transition. 

In  Fig.  4  we  plot  the  same  FFT  images  as  in  Fig.  3  at  the  temperature  T=1 100K.  The  lattice  structure  of  the 
top-most  layer  is  still  maintained  at  this  temperature,  but,  with  a  careful  inspection  (especially,  at  second-order  peaks 
in  the  upper  half  of  the  graph),  the  peaks  from  the  rotated  domains  are  significantly  suppressed.  The  rotated  domains 
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become  unstable  at  high  temperature,  and  an  orientational  phase  transition  from  the  rotated  domain  to  unrotated 
domain  at  temperature  around  T  ~  1000K  is  expected. 

4.  Conclusion.  We  have  studied  thermodynamic  behavior  of  the  hexagonal  top-most  layer  of  Au  (001) 
surface  using  a  MD  simulation  technique.  We  calculate  the  interlayer  distances  as  a  function  of  temperature.  We 
find  that  the  layer  distance,  dj,  between  the  top-most  layer  and  second  layer  is  contracted  as  the  temperature 
increases  whereas  the  second  layer  distance,  d2,  decreases  or  increases  depending  on  the  cluster  size.  The  thermal 
expansion  coefficients  are  —  3.05  x  10-5  A  /  K  for  (101x87)  A  and  —4.93  x  10  5  AIK  for  (262  X227)A  clusters. 
The  FFT  images  of  the  atomic  density  show  that  the  rotated  domains  with  rotational  angles  9=1°  coexist  with  the 
unrotated  domains  in  the  top  layer.  As  temperature  increases  these  rotated  domains  undergo  a  phase  transition  into 
unrotated  domains.  Our  findings  are  in  very  good  agreement  with  experiments  and  previous  theories. 
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